The Centrality Dependence of the Parton Bubble Model for high energy heavy 
ion collisions and fireball surface substructure at RHIC 
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In an earlier paper we developed a QCD inspired theoretical parton bubble model (PBM) for 
RHIC/LHC. The motivation for the PBM was to develop a model which would reasonably quanti- 
tatively agree with the strong charged particle pair correlations observed by the STAR collaboration 
at RHIC in Au + Au central collisions at ^/snn ~ 200 GeV in the transverse momentum range 0.8 
GeV/c to 2.0 GeV/c. The model was constructed to also agree with the Hanbury Brown and Twiss 
(HBT) observed small final state source size ~2fm radii in the transverse momentum range above 
0.8 GeV/c. The model assumed a substructure of a ring of localized adjoining ~2fm radius bubbles 
perpendicular to the collider beam direction, centered on the beam, at mid-rapidity. The bubble ring 
was assumed to be located on the expanding fireball surface of the Au -I- Au collision. These bubbles 
consist of almost entirely of gluons and form gluonic hot spots on the fireball surface. We achieved 
a reasonable quantitative agreement with the results of both of the physically significant Charge 
Independent (CI) and Charge Dependent (CD) correlations that were observed. In this paper we 
extend the model (PBME) to include the changing development of bubbles with centrality from the 
most central region where bubbles are very important to the most peripheral where the bubbles are 
gone. Energy density is found to be related to bubble formation and as centrality decreases the 
maximum energy density and bubbles shift from symmetry around the beam axis to the reaction 
plane region causing a strong correlation of bubble formation with elliptic flow. We find reasonably 
quantitative agreement (within a few percent of the total correlations) with a new precision RHIC 
experiment which extended the centrality region investigated to the range 0-80% (most central to 
most peripheral). The characteristics and behavior of the bubbles imply they represent a significant 
substructure formed on the surface of the fireball at kinetic freezeout. 

PACS numbers: 25.75.Gz, 12.38.Mh 



I. INTRODUCTION AND REVIEW OF 
MODEL 

In the early eighties L. van Hove[l| proposed a 
bubble model as a way of finding convincing evi- 
dence for a Quark-Gluon Plasma (QGP). His model 
based on a string calculation predicted that in one 
to perhaps a few observable localized rapidity bumps 
would appear in the final state of some events and he 
gave a prescription for experimentally finding them. 
We searched for these in all relevant data, but we nor 
anyone else ever found experimental evidence for the 
rapidity bumps. 

There were numerous other bubble models and 
some charge correlation models proposed. Some ex- 
amples of these are references [2, 3,, A, Q, 7] . How- 
ever, to our knowledge no significant experimental 
evidence has been found for any of these. 

In a previous publication Q we proposed a parton 
bubble model (PBM) for central (impact parame- 
ter near zero) high energy heavy ion collisions at 
RHIC/LHC which contains a substructure consist- 
ing of an 8fm radius single ring of a dozen adjoining 
2f radius bubbles transverse to the collider beam di- 
rection, centered on the beam, and located at or near 
mid-rapidity on the expanding surface of the fireball 



at kinetic freezeout. The bubble radius, and the 
bubble ring radius were estimated_considering the 
Hanbury-Brown and Twiss (HBT 9]) observations, 
and other general considerations utilizing the blast 
wave model. We assumed these bubbles (gluonic hot 
spots) are likely the final state result of quark-gluon 
plasma (QGP) formation. Thus this is the geometry 
for the final state kinetic freezeout of the QGP bub- 
bles on the surface of the expanding fireball treated 
in a blast wave model. In the central (near impact 
parameter 0) mid-rapidity region at RHIC we are 
observing the region where the highest energy den- 
sities and temperatures (parton energies) are pro- 
duced. The experimentally observed y/sJ/N = 200 
GeV central Au -I- Au collisions at RHIC|Llfl|] produce 
initial energy densities [ll| which exceed those pre- 
dicted by lattice quantum chromodynamics (QCD) 
as sufficient for production of a QGP[T^. 

This single bubble ring resides at mid-rapidity on 
the surface of the expanding fireball at kinetic freeze- 
out. Thus each bubble would emit a considerable 
fraction of final state particles (observed experimen- 
tally) resulting from the QGP state. The fraction 
of all the final state particles from bubbles is ~ 1/2. 
There would be very little re-interaction for parti- 
cles emitted outward from the surface because the 
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final state surface of the fireball is at kinetic freeze- 
out. The bubble substructure (surrounded by cooler 
background) results in the lumpy surface of the fire- 
ball at kinetic freezeout. Section II presents more 
detail on the assumptions made, the development, 
and construction of the PBM. 

The PBM successfully explained in a reasonably 
quantitative manner all of the particle pairs cor- 
relations in a precision STAR central production 
experiment [To|. See Section 4 of Ref-.S] for this 
comparison. In Ref.Q some aspects of quark- 
quark recombination were compared with the PBM 
and a good agreement was obtained (see Section 5 
PBMQ). 

In this paper we extend the model of Ref. Q to 
consider the case of varying the centrality bins (Sec- 
tion III). We wish to compare with a new RHIC 
y^sjvjv = 200 GeV minimum bias Au -I- Au analysis 
which covered the 0-80% centrality range p^. 

This paper is organized as follows: 

Section I is the Introduction. Section II sum- 
marizes the assumptions made in the prior model 
(PBM) for the central region (0-10% centrality^]) 
Section II also discusses the relevance and the rea- 
soning behind these assumptions. Section III dis- 
cusses extending the PBM to the Parton Bubble 
Model Extended (PBME) so that it becomes able 
to reasonably quantitatively fit and explain the new 
0-80% centrality high precision data[l3|. Section IV 
discusses general characteristics of the PBM. Section 
V presents and discusses a comparison of the ex- 
perimentally determined Charge Independent (CI) 
correlation with PBME as a function of % central- 
ity. Section VI presents and discusses a comparison 
of the experimentally determined Charge Dependent 
(CD) correlation with PBME as a function of % cen- 
trality. Section VII presents further details of the 
PBME bubble correlation. Section VIII is the Sum- 
mary and Discussion. 



II. ASSUMPTIONS AND DEVELOPMENT 
OF THE PBM 

Our goal for the past 2 decades was to develop 
a model of bubble production in relativistic heavy 
ion collisions assumed to originate from a QGP, 
which could be reasonably quantitatively compared 
with relevant experimental data of sufficient preci- 
sion and scope. Thus hopefully we could obtain con- 
vincing, or at least substantial evidence for the ex- 
istence of bubble production. Since this is a pro- 
cess which clearly involves strong non-perturbation 
QCD, pQCD calculations can serve only as a rough 
guide. 

Thus we concluded we needed to obtain a strong 



hint from experimental data as evidence for possi- 
ble bubble substructure to proceed to build a re- 
alistic model. The failures of obtaining significant 
experimental evidence for the many bubble models 
which did not incorporate such a strong hint from 
the experimental data led us to conclude that it was 
essential to obtain one from experimental observa- 
tions. Then one could build the bubble model based 
on observations in order to provide quantitative de- 
scription of the data. 

We utilize a two particle correlation function in 
the two dimensional (2-D) space of vs Ary. The 
azimuthal angle (/) of a particle is defined by the angle 
of the particle with respect to the vertical axis which 
is perpendicular to the beam axis and is measured 
in a clock-wise direction about the beam. A(/) is the 
difference, 0i - 02 , of the (j) angle of a pair of particles 
(1 and 2). The pseudo-rapidity 77 of a particle is 
measured along one of the beam directions. A77 is 
the difference, rji - 772, of the ry values of a pair of 
particles (1 and 2). 

The two dimensional (2-D) total correlation func- 
tion is defined as: 



C(A0, Ary) = 5(A0, A77)/M(A0, A77). (1) 

Where S(A0, At?) is the number of pairs at the 
corresponding values of A0, At; coming from the 
same event, after we have summed over all the 
events. M(A0, At;) is the number of pairs at the cor- 
responding values of A^, A77 coming from the mixed 
events, after we have summed over all our created 
mixed events. A mixed event pair has each of the 
two particles chosen from a different event. We make 
on the order of ten times the number of mixed events 
as real events. We rescale the number of pairs in the 
mixed events to be equal to the number of pairs in 
the real events. 

The behavior of the HBT quantum interference 
radii 0, especially Rs was interpreted as indicative 
of observation of spatial radii ~2fm in RHIC central 
Au -I- Au collisions at ^snn ~ 200 GeV for pt > 
0.8 GeV/c. Rs reduced to ^2fm from ^6fm as pt 
increased from ~ 0.2 GeV/c to greater than ~ 0.8 
GeV/c. HBT quantum interference is by necessity 
only measured for pairs of charged particles of the 
same sign with small difference in momentum. In 
our model all particles coming from a bubble are 
constrained to come from a ~ 2fm radius. Particles 
above 0.8 GeV/cpt that come from different bubbles 
differ in momentum such that they do not show HBT 
quantum interference. 

The generally accepted explanation for observing 
these increasingly smaller final state HBT radii as 
Pt increased was that radial flow increasingly fo- 
cused the viewed region of the overall final state 
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source into a smaller volume as pt increased. We 
refer to this phenomenon as phase space focusing 
due to flow. We realized (in our first bubble model 
paper [U published in 2003) that flow phase space 
focusing implies that the viewed region on the sur- 
face for > 0.8 GeV/c of the fire ball would have 
a volume with radii of ~2fm. With in this volume 
a hot source producing a larger number of particles 
would move out away from the surface, and being 
focused together would lead to an increase of parti- 
cles emitted in an angular region. However the HBT 
correlation function has the property that a ring of 
essentially similar bubbles as assumed in our PBM 
model would image on top of each other forming an 
average bubble which HBT would be viewing. The 
lower Pt cut of 0.8 GeV/c would allow HBT to view 
and resolve this average bubble formed from differ- 
ences in momentum of two particles for pt above 0.8 
GeV/c. Our model populated these HBT viewed 
regions with a ring of bubbles or gluonic hot spots 
producing a larger number of particles with an angu- 
lar correlation. We form a correlation function based 
on the difference of angles between two charged par- 
ticles which images the 12 bubbles on top of each 
other. 

It is important to note that we have calculated 
that the correlations observed at RHIC are strong 
enough so that if there was only one bubble instead 
of a ring of bubbles distributed around the azimuth 
as our PBM assumed, one would produce an angu- 
lar region with huge amplitude spikes in individual 
events. These spikes are not observed in individ- 
ual events at RHIC. Therefore the RHIC correlation 
data has to be built of smaller distributed correlated 
regions as assumed by the PBM. 

We had used a virtually identical bubble ring and 
the same correlation function in our first bubble 
modeim which was published in 2003 and was able 
to predict and subsequently explain important gen- 
eral characteristics of the experimental analysis of 
the two charged particle correlation data for central 
collisions [l0|. This experimental paper used the pt 
cuts developed in our 2003 paper 0.8 < pt < 2.0 
GeV/c. The 2 GeV/c cut was employed to make 
jet contamination negligible. The angular correla- 
tions used the differences of the azimuthal angles 
(A(/)) and the differences of psucdorapidity (Ajy) of 
all charge particle pairs which had the property that 
imaged all bubbles in the ring on top of each other. 
This allowed the phase space focusing by radial flow 
to provide consistency with the observed correlation 
data analyses and the HBT observations. 

An important assumption made in our PBM is 
that in a central heavy ion collision (e.g. Au -I- 
Au) at high RHIC energies a high density of en- 
ergetic partons (virtually all gluons) form a dense 



opaque fireball. This dense opaque fireball has a 
large amount of radial flow and can be described very 
well by a blast wave model. The usually employed 
blast wave model we usedQ has its maximum veloc- 
ity at the surface of the fireball at kinetic freezeout 
of approximately (3/4)c. 

A theoretical pQCD calulationfTsj in 1987 con- 
cluded that jets formed with initial parton transverse 
momenta of around 3 GeV/c (also applicable down 
to 2 GeV/c) would become thermalized in a ^/sWn 
= 200 GeV \J + 1] collision at RHIC and would not 
escape from the system. Therefore these jets would 
not result in the correlations observed at RHIC . 
One could speculate that these thermalized jets form 
the dense opaque fireball. 

There is direct experimental evidence of strong 
quenching of high pt particles at RHIC (e.g. 
Refs.flGl. IitI. [l^). An experimental result which is 
independent of our model demonstrates that in the 
central region charged particles with pt of (0.8 < pt 
< 4.0 GeV/c) are emitted from the surface of the 
dense opaque fireball. This point was demonstrated 
in Section IV B of Ref.[l3l. We quote the last 3 
sentences from this Section IV B which states the 
conclusion. "This surface or near surface hadroniza- 
tion and emission from the fireball both occur in the 
central region and all other centralities where there 
is appreciable particle density. In the most periph- 
eral bins the particle density is low enough to allow 
undisturbed fragmentation and thus no change in 
the CD correlation. Thus the CD behavior is consis- 
tent with a surface emission model such as Ref.Q." 
This result implies that surface or near surface emis- 
sion occurs for charged particles with pt (0.8 < pt < 
4.0 GeV/c) for all centralities. 

The HIJING event generator 19j combines Pythia 
jetsjl^l and the Lund model[21| and thus was a fa- 
miliar base for constructing our PBM model. As 
we had done in our earlier bubble model(l3|. we re- 
placed the Pythia jets in HIJING with our bubble 
ring. Momentum, energy and charge conservation 
are all satisfied within the bubble ring in the PBM. 
The bubble ring becomes the source of cmmitted 
particle correlation generated by Pythia which we 
used for fragmentation of the bubbles. We made 
the approximation that hard jet particles are essen- 
tially removed or have their correlations removed by 
quenching. The only remaining particles from HI- 
JING are the beam jet fragmentation particles which 
are soft and have no correlation. These particles be- 
come our background particles in the PBM. How- 
ever we did include the effects of elliptic flowfl^l on 
the soft beam jet fragmentation particles, since el- 
liptic flow does generate a small cos(2A(/>) term in 
the correlation even in the central collisions. The 
procedure we employed was to in each event deter- 
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mine the reaction plane and modulate the soft beam 
fragmentation particles by the elliptic flow term 2v2 
cos(2A(/)) which was a sufficient approximation for 
elliptic flow effects. Thus the elliptic flow effects 
were put into the model on an event by event basis. 
In the PBM we used angular correlations of charge 
particle pairs to predict and fit the experimentally 
determined correlations [lOj . 

In Fig. 3 of the PBM publication [1] we schemati- 
cally show a number (3-4) of parallel parton showers 
contained in each bubble. The parton showers must 
be parallel to give results which are consistent with 
the experimental analysis of angular correlations. 

The particle production from our bubbles uses a 
similar parton QCD shower fragmentation as a jet 
with a well defined <j) angle (Fig.2 of the PBMQ). 
The pt distribution of the partons inside the bub- 
ble is similar to pQCD but has a suppression in 
the high pt region like the dataliet [ij, The 
3-4 partons have different longitudinal momentum. 
At kinetic freezeout we used Pythia fragmentation 
functions [20] for the bubble fragmentation to form 
the final state emitted charged particles. 

Models that successfully predict and fit non- 
perturbative QCD experimental results reasonably, 
almost always have to adjust some parameters when 
comparing with experimental data. In the PBM 
there are two such adjustable parameters which are 
the number of partons in a bubble and the longitudi- 
nal momenta of the partons. Ever since Landau dis- 
covered in cosmic rays a long time ago that excited 
nuclear fireballs exhibited a longitudinal expansion 
this fact was well known. Adjusting the longitudi- 
nal momenta of the partons is obviously necessary 
to explain the expansion in A77 in the central pro- 
duction experiment . The two parameters in the 
PBM were adjusted for a set number of bubbles by 
comparing the PBM fit to the experimentally deter- 
mined final state charge pair correlation (CI defined 
below) . 

In our models (e.g. PBM) the experimental cor- 
relation data analysis we are comparing to utilizes 
the correlations of charge particle pairs. There are 
2 types of such correlated pairs, namely unlike-sign 
charge pairs (US) and like -sign charge pairs (LS). 
The total sample of correlated charge pairs gener- 
ated and emitted in the final state of a theoretical 
model is equal to the average of the US and LS cor- 
relations ((US -I- LS)/2). In an experiment the total 
sample detected depends on the acceptance and ef- 
ficiency of the detector. One makes cuts on the the- 
oretical model to account for these effects. However 
the correlation function used in equation 1 is conven- 
tionally used in experimental analyses and is drasti- 
cally independent of acceptance effects as stated in 
Section 2.l\§\ allowing reasonable comparisons. 



The Charge Independent (CI) correlation is con- 
ventionally defined as the unlike-sign charge pair 
correlation (US) -I- like-sign charge pair correlation 
(LS). The total correlation derived when using all 
particle pairs independent of what charge signs are 
used to form the correlation is equal to CI/2. Thus 
CI/2 gives the average structure of the correlated 
emitting sources independent of charge and repre- 
sents the overall physical phenomenon. 

The Charge Dependent (CD) correlation is con- 
ventionally defined as the US - LS. The subtrac- 
tion of the total like-sign charge pairs correlation 
(LS) in forming the CD is equivalent to removal of 
the opposite sign charge pairs which are not from 
the same space-time region where charge has to 
balance [Tol H^. Therefore the CD is expected to 
represent the correlation of unlike-sign pairs from 
the same space-time region where charge is balanced 
as modified by interaction with the medium. In our 
surface emission model (PBM) the interaction with 
the medium is absent therefore, the CD is expected 
to exhibit the correlation of Pythia jets which are 
produced in a vacuum. 

In Section 4.2 of the PBMQ we compare the total 
CI correlations of PBM and precision experimental 
correlation data[l3]. We adjusted the two param- 
eters by comparing the CI of the PBM with the 
CI of the data. We show the PBM and the data 
agree within less than 10% of the total CI correla- 
tion which is a reasonably quantitative agreement. 

In Section 4.3 of the PBM 8] we compare the CD 
correlations of PBM with the experiment by making 
use of the fact that the net charge fiuctuation sup- 
pression is directly related to an integral over the 
CD. Thus we compared the net charge fluctuation 
suppressions of the data with the PBM and found 
agreement within errors. We did the comparison of 
the CD in this way since the experimental paper we 
compared the model with considered this to be an 
important aspect of the CD correlation and chose to 
treat the CD in this manner. 



III. EXTENSION OF PBM TO COVER 
0-80% CENTRALITIES 

The PBM is a successful bubble model for central 
heavy ion (e.g. Au -I- Au) collisions at the highest 
energy at RHIC as discussed in the previous section. 
However an interesting question that arises is what 
would happen to bubble production and the general 
characteristics of charge pair correlations as central- 
ity varies from most central (impact parameter near 
zero) to peripheral collisions. The PBM has been 
successfully tested in the centrality range approxi- 
mately (0-10%). A new precision RHIC minimum 
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bias trigger data analysis for Au + Au collisions at 
i/snn = 200 GeV in transverse momentum range 
0.8 GeV/c to 4.0 GeV/c covers the 0-80% centrality 
range flSj and is ideally suited for investigating the 
varying centrality evolution of the PBM. This exper- 
imental analysis was done in a manner which was a 
logical extension of the central production paper |lC|. 
The data, the charge particle pairs correlations (US 
and LS), the CI and the CD were treated in a similar 
manner, but analyzed and fit separately in each of 
nine centrality bins. 

In the data analysis p!^ the most central bins US, 
LS, and CI were consistent with the results of the 
prior central production experiment. However as 
one moves from central to peripheral bins a jet-like 
component is increasingly evident (in the data anal- 
ysis) till the most peripheral bins where there is only 
jet-like correlations. The elliptic flow amplitude 2?;| 
as part of the correlation increases as one moves 
to more peripheral centralities. This flow reaches 
its maximum at 40-50% centrality (see Fig. 28 of 
Ref.dl). 

We will proceed to extend the PBM to include 
the entire 0-80% centrality region. We name this 
extended version of the PBM as PBME. Several ob- 
vious characteristics of the experimental and theo- 
retical analysis that must be included in this exten- 
sion of the model to PBME are discussed below. 

The previous central collision bubble ring geome- 
try was well suited for the most central collisions sit- 
uation since the highest energy densities are circular 
around the beam which is the expected geometric 
symmetry. As one moves from central toward pe- 
ripheral the decreasing size and change of shape of 
the overlap region of the two Au nuclei determines 
where the energy densities are highest. The overlap 
of matter in the two Au nuclei becomes greater in 
the reaction plane region, while the overlap of mat- 
ter becomes less outside of the reaction plane. This 
breaking of symmetry will modify the overall spatial 
shape and location of the bubbles. 

The effects of elliptic flow '22'] on the events in each 
centrality bin were put into the model using the same 
procedure we used as previously described in Section 
II for the PBM. The procedure we employed was to 
in each event determine the reaction plane and mod- 
ulate the soft beam fragmentation particles by the 
elliptic flow term 2v2 cos(2A0) which was a suffi- 
cient approximation for elliptic flow effects. Thus 
the elliptic flow effects were put into the model on 
an event by event basis. 

Jet quenching is largest in the most central colli- 
sions and decreases as one moves to more peripheral 
bins. We found a sufRcient way to put the effects of 
strong jet quenching in the central collisions. We set 
jet quenching to its maximum in HIJING for 0-30% 



centralities by removing all jets (quenched away). 
For centralities 30-80% we use the non-jet-quenching 
version of HIJING thus all jets become part of the 
event. The soft beam jet particles have elliptic flow 
as described above. 

We relied on the blast wave to determine the 
geometry of where the energy density is highest. 
These regions of high energy density is where bub- 
bles are formed. In the final state at kinetic freezeout 
the bubbles are located on the fireball (blast wave) 
surface which is the source of emmitted correlated 
charge particle pairs generated by Pythia fragmenta- 
tion of the bubbles. Thus the blast wave determined 
the location, number and geometry of the bubbles. 



IV. GENERAL CHARACTERISTICS OF 
THE PBME 

For the most central bin we have the same bubble 
geometry and partons per bubble as Ref.[^]; namely 
12 bubbles and 3-4 partons per bubble. As we move 
away from the most central the number of partons 
per bubble decreases dropping to 3 then 2 and fi- 
nally 1 at 50-60% centrality. The number of bubbles 
formed per event decreases from 12 in the most cen- 
tral bin to 0.3 in the 50-60% bin (see Table I in this 
section). 

The most central collisions have a bubble ring 
symmetry about the beam axis (see Fig.l of Ref. Q). 
The region of highest energy density will also be 
symmetric about this axis and the ring of bubbles 
is the expected geometry. As we move to more pe- 
ripheral collisions the symmetry becomes defined by 
the reaction plane. The region of highest energy den- 
sity becomes more concentrated in the region of the 
reaction plane. Therefore it is reasonable to expect 
the ring symmetry will be broken and the bubbles 
farthest from the reaction plane would disappear. 
We find comparing 0-5% centrality to 5-10% cen- 
trality the 12 bubbles reduce to 10 bubbles with the 
bubbles perpendicular to the reaction plane gone. 
For 10-20% the number drops to 7 bubbles with 4 
bubbles being near the reaction plane and 3 appear- 
ing above and below the plane, while for 20-30% we 
drop to 5 bubbles with 4 bubbles being near the re- 
action plane with only 1 appearing above or below 
the plane. Moving to 30-40%, 40-50%, and 50-60% 
the bubbles are in the reaction plane region with the 
probability per event of making a bubble being 1.5, 
0.6, 0.3 respectively (see Table I). The above changes 
in bubble production are due to the decrease in en- 
ergy density as the % centrality range becomes more 
peripheral which decreases the overlap region of the 
Au -I- Au colliding nuclei. 

The STAR experiment has measured charged par- 
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tide pair correlations for minimum bias An + Au 
events at ^/SNN — 200 GeVjJI]- The pt range of 
that data is 0.8 to 4.0 GeV/c for the entire - 360° 
(j) range and the 77 range |r/| < 1.0. Both the exper- 
imental data analysis compared to and the PMBE 
model utilize two dimensional A(f) Arj correlations. 
See examples of these two dimensional perspective 
plots in Refs.[i,[l3| and Fig.11-12 (this paper). 

In order to compare and present these two dimen- 
sional plots for the new data[3] and the PBME we 
divided the entire A77 region into five Arj bins which 
covered the entire Arj range. Each Arj bin could then 
be presented as a one dimensional projection and a 
comparison between data and model can be made. 
The five Arj bins were 0.0 to 0.3, 0.3 to 0.6, 0.6 to 0.9, 
0.9 to 1.2 and 1.2 to 1.5. In each of these five Aij bins 
the A(j) correlations for the charged particles covered 
the entire Acj) range - 180°. Due to demonstrated 
symmetry in the data and the model the 360° A(p 
range which was experimentally detected was folded 
resulting in the 180° ranges. 



Table I 


Centrality 


bubbles per event 


partons per bubble 


0-5% 


12 


3-4 


5 - 10% 


10 


3 


10 - 20% 


7 


3 


20 - 30% 


5 


3 


30 - 40% 


1.5 


2 


40 - 50% 


0.6 


2 


50 - 60% 


0.3 


1 



Table I. Parameters of bubble model with cen- 
trality. 



V. CHARGE INDEPENDENT (CI) 
CORRELATION 

In this section we are going to compare the PBME 
model predicted CI correlations with the STAR ex- 
perimental analysis results [13]. The two particle 
correlations are formed from two different types of 
charge particle pairs: 

1) Unlike-Sign charge pairs (US) 

2) Like-Sign charge pairs (LS) 

The Charge Independent (CI) total correlation = 
US -I- LS. It is the sum of the CI signal -I- Back- 
ground correlations in the final state after kinetic 
freezeout. Thus it is the total two particle correla- 
tion observed in the detector (STAR TPC). The en- 
tire CI correlation (signal -I- background) is used for 
comparing the analysis results with the model. This 
eliminates any model dependence on the separation 



of signal from background. Certain necessary cor- 
rections and cuts in the experimental analysis were 
applied to the PBME model (this paper) so that 
a quantitative comparison could be made with the 
experimental analysis. The CI displays the average 
correlation structure of the emitting sources in the 
final state after kinetic freezeout and thus is physi- 
cally significant. 

The five plots comparing the experimental analy- 
sis of the CI and the PBME fits are shown in Fig. 
1-5. On the vertical axis the CI is multiplied by the 
average event multiplicity within the particular cen- 
trality bin shown by the symbols on the plot. This 
procedure is necessary when one compares different 
centralities in order to make the comparison inde- 
pendent of multiplicity. 

The correlation function is given in equation 1. 
The experimentally observed correlation has a nu- 
merator which is proportional to the number of cor- 
related particle pairs which itself is proportional to 
the multiplicity. However the denominator is pro- 
portional to the total number of pairs which can be 
formed which is proportional to the square of the 
multiplicity. Therefore in order to make compar- 
isons of different centralities or other experiments 
one multiplies by the average multiplicity to remove 
the dependence on the multiplicity. This procedure 
is referred to as multiplicity scaling or multiplicity 
scaled. 

For each centrality we generated 500,000 simu- 
lated events with the bubble geometries presented 
above. The impact parameter range of HIJING for 
the different centralities is given in Table II. See Ta- 
ble I for the number of bubbles per event and partons 
per bubble as a function of % centrality bins. 



Tabic II 


Centrality 


Impact parameter (fm) 


0-5% 


0.0-2.9 


5 - 10% 


2.9-4.1 


10 - 20% 


4.1-5.8 


20 - 30% 


5.8-7.1 


30 - 40% 


7.1-8.2 


40 - 50% 


8.2-9.2 


50 - 60% 


9.2-10.1 


60 - 70% 


10.1-10.9 


70 - 80% 


10.9-11.7 



Table II. Impact Parameters of HIJING with 
Centrality. 

In Fig.l to Fig. 5 we used the STAR calculated 
Charge Independent (CI) correlation which is the 
average correlation of the unlike-sign charge pairs 
correlation plus the like-sign charge pairs correla- 
tion. In order to allow comparison with the different 
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centrality bins we and also STAR multiplied by the 
multiplicity, because as explained previously this re- 
moves the dilution of the signals due to the quadratic 
increase of pair combinations. The multiplicity that 
is used in Ref.fTs'] is based on the particles mea- 
sured in the STAR TPC (Time Projection Cham- 
ber). There are readout bomidaries between the 12 
TPC sectors which cover the azimuth that do not 
measure tracks. We put these readout boundaries in 
our Monte Carlo generation of our particles, which 
cause a loss of approximately 10% of the particles. 

The multiplicity scaled correlation for each cen- 
trality for a given Ar? bin is plotted with the maxi- 
mum angle of the away side (Ac/) = 180°) shifted to 
the same value. The horizontal line for each central- 
ity shows the shifted average multiplicity line which 
was normalized to a mean of 1 in the original CI 
correlation and became equal to the average mul- 
tiplicity after becoming rescaled to the scaled cor- 
relation. The solid curves are the PBME calcula- 
tions shifted by the same amount as the data. The 
agreement between the PBME and the RHIC data in 
Fig.l to Fig. 5 is within a few percent of the total cor- 
relation. Considering that the model (PBME) does 
not completely include important non-perturbative 
QCD effects contained in the data we consider this 
a reasonable quantitative agreement. The CI x mul- 
tiplicity displays the average structure of the corre- 
lation sources at kinetic freezeout. 



VI. CHARGE DEPENDENT (CD) 
CORRELATION 

We compare total experimentally observed CD 
correlations to total theoretically predicted CD cor- 
relations in order to avoid possible uncertainties due 
to separation of signals and background. 

The Charge Dependent (CD) correlation which is 
the difference between the unlike-sign charge pair 
correlations and the like-sign charge pair correlations 
( US - LS ) displays a measure of the emission cor- 
relation of the opposite sign pairs of particles emit- 
ted from the sam e sp ace-time region at the time of 
hadronizationfiol l23j| . The CD as all correlations in 
all centrality ranges in the PBME (which includes 
a HIJING jet component) has it's fragmentation of 
all partons determined by Pythia fragmentation fl^. 
We use the same projection method of the 2-D CD 
correlation into the same five A77 ranges and multi- 
ply the CD by the multiplicity as described and dis- 
cussed in the previous Section V. In Fig. 6 we show 
the multiplicity times the CD A(/) correlation from 
Ref.[l3| within the A?/ range 0.0 to 0.3 for each cen- 
trality compared to the bubble model (PBME) cal- 
culations. We achieve a good agreement between 



data and model. Similar results and good agree- 
ment with the PBME occurs for the four other A77 
ranges not shown. Thus we have good agreement in 
all 5 Ary ranges which together comprise the entire 
CD. 

In Fig.7 we plot from Ref.fTs"] the CD A0 corre- 
lation for 4 A77 bins covering the range 0.0 < At] < 
1.2.^ Each of 9 centralities shown were scaled so 
that the 5 - 10° A^ bin for 0.0 < Ai] < 0.3 range is 
normalized to 1. This was done in order to remove 
the scale difference between the different centrality 
ranges and allow us to show that the CD shape is ap- 
proximately independent of centrality. In Fig.7 the 
experimental analysis points in each A77 range clus- 
ter around the four lines which corresponds to each 
of the four Ary ranges generated by Pythia iets[20j. 
Pythia jets were used as the jets in HIJING and 
Pythia fragmentation was used in the parton bubble 
model @. Thus the CD shape is in good agreement 
with Pythia and is independent of centrality. Pythia 
jet CD correlations are the initial correlations of op- 
posite sign charge pairs from the same space-time 
region at the time of hadronization[l3,[2^. The fact 
that there is essentially no change in these correla- 
tions at the time of kinetic freezeout demonstrates 
that there is little or no further interaction of these 
opposite sign charge pairs with the fireball medium 
from the time of hadronization till the time of ki- 
netic freezeout. Thus both hadronization and kinetic 
freezeout occur at or very near the surface of the ex- 
panding fireball at all centralities. Hence in general 
the fireball is dense and opaque at kinetic freezeout. 
One should note that in the peripheral bins there 
is very little matter in the short path length from 
any point to the surface. This is consistent with sur- 
face emission from the fireball at kinetic freezeout or 
undisturbed fragmentation in those peripheral bins 
where the path length to the surface is small. 

Both the original parton bubble model (PBM)0] 
and the present extension of the model (PBME) con- 
struct surface bubbles which are boosted by the ex- 
panding fireball. These bubbles at freezeout give 
results that are consistent with experimental corre- 
lation data[ll|- Furthermore the PBME fits the ob- 
served CD correlation to within a few percent of the 
correlation at all centralities. The PBME produces 
the CD shape that is consistent with Pythia jet CD 
correlations. As previously shown in Section V the 
PBME fits the observed CI correlations to within 
a few percent of the correlation. The US and LS 



^ The Ari range 1.2 to 1.5 (not shown) is just a flat back- 
ground with very Uttle CD signal from Pythia parton frag- 
mentation left. 
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o70 % to 80 % 

o,n I I \ I I \ I I \ I I \ I \ \ I \ \ 

20 40 60 80 100 120 140 160 180 

A(|) (deg) 

FIG. 1: The multiplicity(MULT) times the CI correlation vs. for 0.0 < Ar; <0.3. Nine centralities are shown from 
70% to 80% with centrality increasing to 0% to 5%. 216 is the multiplicity for 0% to 5% and all other centralities 
are shifted up so that the 180° value is equal. Each multiplicity for each centrality is shown shifted. The solid curves 
are the bubble model (PMBE) calculations shifted by same amount as the data. The PBME and RHIC data agree 
within a few percent of the total CI correlation in each centrality bin in each Arj range for all A<j} angles. 



correlations are linear combinations of the CI and 
CD. Therefore they arc also fit to within a few per- 
cent of the correlation. Thus all two charge particle 
pair correlations are reasonably quantitatively fit by 
the PBME (Parton Bubble Model Extended). One 
should note that the fireball in both the PBM and 
PBME is treated in a blast wave model with the 
bubbles forming on the surface and emitting their 
final state particles at kinetic freezeout. The fireball 
surface is moving at the maximum velocity at kinetic 
freezeout. 



VII. FURTHER DISCUSSION OF THE 
PARTON BUBBLE MODEL CORRELATION 



We note from Tables 1-111 that as we move from 
the most central region (impact parameter small) 
toward the peripheral bins the number of bubbles in 
a centrality bin per event and the number of partons 
per bubble both tend to decrease. In the two most 
peripheral bins bubble production is negligible. 
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o70 % to 80 % 

o,n I I \ I I \ I I \ I I \ I \ \ I \ \ 

20 40 60 80 100 120 140 160 180 

A(|) (deg) 

FIG. 2: The multiplicity(MULT) times the CI correlation vs. for 0.3 < Ar/ <0.6. Nine centralities are shown from 
70% to 80% with centrahty increasing to 0% to 5%. 216 is the multiplicity for 0% to 5% and all other centralities are 
shifted up so that the 180° value is equal. Each multiplicity for each centrality is shown shifted. The solid curves are 
the bubble model (PBME) calculations shifted by same amount as the data. The PBME and the RHIC data agree 
within a few percent for the total CI correlation 



Table III 


No partons 


particles per bubble 


Pt per bubble(GeV/c) 


energy per bubble (GeV) 


3-4 


8.5 


7.0 


9.8 


3 


7.0 


5.7 


8.0 


2 


4.9 


4.0 


4.9 


1 


2.6 


2.1 


2.1 



Table III. Parameters of bubble model per num- 
ber of partons. 
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MULTxCI for 0.6 < Ar| < 0.9 




0.0% to 5.0% 
•5.0% to 10 % 

aIO % to 20 % 
t20 % to 30 % 
o30 % to 40 % 
□40 % to 50 % 
a50 % to 60 % 
o60 % to 70 % 
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60 



80 



100 



120 



o70 % to 80 % 

J \ I \ 

140 160 



180 



A(|) (deg) 



FIG. 3: The multiplicity(MULT) times the CI correlation vs. for 0.6 < Ar/ <0.9. Nine centralities are shown from 
70% to 80% with centrality increasing to 0% to 5%. 216 is the multiplicity for 0% to 5% and all other centralities are 
shifted up so that the 180° value is equal. Each multiplicity for each centrality is shown shifted. The solid curves are 
the bubble model (PBME) calculations shifted by same amount as the data. The PBME and the RHIC data agree 
within a few percent for the total CI correlation 
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MULTxCI for 0.9 < Ar| < 1.2 




0.0% to 5.0% 
•5.0% to 10 % 

aIO % to 20 % 
t20 % to 30 % 
o30 % to 40 % 
□40 % to 50 % 
a50 % to 60 % 
o60 % to 70 % 
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140 160 



180 
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FIG. 4: The multiplicity(MULT) times the CI correlation vs. Acf} for 0.9 < Ar/ <1.2. Nine centralities are shown from 
70% to 80% with centrality increasing to 0% to 5%. 216 is the multiplicity for 0% to 5% and all other centralities are 
shifted up so that the 180° value is equal. Each multiplicity for each centrality is shown shifted. The solid curves are 
the bubble model (PBME) calculations shifted by same amount as the data. The PBME and the RHIC data agree 
within a few percent for the total CI correlation 
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FIG. 5: The multiplicity(MULT) times the CI correlation vs. for 1.2 < Ar/ <1.5. Nine centralities are shown from 
70% to 80% with centrality increasing to 0% to 5%. 216 is the multiplicity for 0% to 5% and all other centralities are 
shifted up so that the 180° value is equal. Each multiplicity for each centrality is shown shifted. The solid curves are 
the bubble model (PBME) calculations shifted by same amount as the data. The PBME and the RHIC data agree 
within a few percent for the total CI correlation 
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(deg) 



FIG. 6: The multiplicity(MULT) times the CD correlation vs. for 0.0 < A?; <0.3. Nine centralities are shown 
from 70% to 80% with increasing centrality to 0% to 5%. The solid curves are the bubble model (PBME) calculations 
scaled by the multiplicity. We get good agreement. 



14 



1.3 



1.1 



MULTxCD 
1 



Centrality 




Ar| Range 


■ 0.0% to 5.0 


% 


Range 1 


•5.0% to 10 


% 


0.0 to 0.3 


A 10 % to 20 


% 


Range 2 


▼20% to 30 


% 


0.3 to 0.6 


o30%to 40 


% 


Range 3 


□ 40% to 50 


% 


0.6 to 0.9 


A50%to 60 


% 


Range 4 


60 % to 70 


% 


0.9 to 1.2 


*70%to 80 


% 






80 100 

A(|) (deg) 

FIG. 7: The product of the multiplicity(MULT) and the CD correlation vs. A(j> for 4 Arj bins. Nine centralities are 
shown from 70% to 80% with increasing centrality to 0% to 5%. Each of the centralities were scaled so that the 5 
- 10° A(j> bin for the A-q range 0.0 to 0.3 is normalized to 1. The data points for each Ai] range cluster around the 
Pythia jet predictions (lines). Thus we see that the CD shape is approximately independent of centrality and the CD 
is approximately consistent with Pythia jets[2(?l (line) at each centrality. The PBME fits the observed CD correlation 
to within a few percent of the correlation at all centralities since its predicted shape is given by Pythia jets and we 
get a good agreement for multiplicity scaled CD correlation for each of the centralities. 
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Fig. 8 and Fig. 9 are 2 dimensional perspective 
plots -A^Ayy- of the part of the CI correlation multi- 
plied by the multiplicity (multiplicity scaled) which 
results from particles emmitted from the same bub- 
ble in the centrality bins shown. The A(/) angular 
distributions of the correlation are primarily con- 
fined to the near side (A(/) < 90°). The shape of 
the A0 correlation for 0-30% centrality bins peak 
near small angles and decrease as Acjj increases to- 
ward 90° and are qualitatively similar. However the 
Ary distribution shape remains qualitatively broad in 
the 0-30% centrality bins exhibiting the Atj elonga- 
tion observed in the data[l3. One should note that 
there is very little away side (A0 > 90°) in any of 
these centrality bins. 

Fig. 10 shows a comparison of the bubble signal 
(particle from the same bubble) shape in Arj of the 
different centralities in which we have bubbles. The 
bubbles with the same number of partons are consis- 
tent with the same shape in A77. Thus the longitudi- 
nal (i.e. A77 distribution) determines the number of 
partons in each bubble which determines the longi- 
tudinal shape of the bubbles. We can see that there 
is a maximum expansion in the length of A77 in the 
most central bin (0-5%) and close to this maximum 
is maintained through the (20-30%) bin. However 
the length of Ar/ considerably decreases in the (30- 
60%) bins. 

In Ref. Q the ring of bubbles played an important 
role in the away side (A0 180°) correlation. Corre- 
lation between particles from different bubbles show 
an away side correlation. Fig. 11 and Fig. 12 show the 
correlation that resulted from all particles from all 
the bubbles. We plot this correlation for all central- 
ities where there are bubbles present. 

If we compare Fig. 8 and Fig. 9 (multiplicity scaled 
CI for particles from the same bubble) with Fig. 11 
and Fig. 12 (multiplicity scaled CI for particles from 
all bubbles) we note a striking difference in the away 
side (A0 > 90°) behavior of the two sets. The corre- 
lation resulting from particles emitted by the same 
bubble produces an away side that is very small 
(Fig. 8 and Fig.9). The away side correlation pro- 
duced by particles emitted by all bubbles is large 
in the centrality range 0-30% where bubble pro- 
duction is large and then essentially disappears for 
the more peripheral centralities where the bubble 
production becomes small to negligible (Fig. 11 and 
Fig. 12). This is due to the model conserving mo- 
mentum between the bubbles and is consistent with 
the experimental observations. 

Our model results in bubbles forming where en- 
ergy density is highest. Thus as we have pointed 
out the most central collisions have the highest en- 
ergy density and therefore a circular ring of bub- 
bles perpendicular to and centered on the beam 



axis was able to explain the correlations observed 
in a y/s^N = 200 GeV Au + Au RHIC cen- 
tral production experimental analysis [1, [l0|- The 
FEME model with the same bubble geometry has 
explained the new experimental analysis as a func- 
tion of centrality 13]. However as we move away 
from central production bins toward peripheral bins 
the highest energy density concentrates in the re- 
gion of the reaction plane, the central ring symme- 
try is broken, and the bubbles tend to be produced 
in the reaction plane region. Thus the geometry of 
the bubbles become coupled to elliptic flow. 



VIII. SUMMARY AND DISCUSSION 

In this article we summarize the assumptions 
made, and the reasoning that led to the development 
and construction of our Parton Bubble Modelj^ 
which successfully explained the charge particle pair 
correlations in the central (approximately 0-10% 
centrahty) y/s^N = 200 GeV Au + Au datafioj. 
The PBM was also consistent with the central col- 
lision Au -I- Au HBT results. This is presented and 
discussed in Section 4 of Ref.[^ and Section II of 
this paper. Most of this paper is concerned with ex- 
tending our model which was a central region model 
to be able to treat the geometry of bubble produc- 
tion for 0-80% collision centralities (PBME) such as 
measured and analyzed in recent RHIC data[ll]. In 
the PBME we included elliptic flow for all centrali- 
ties and a jet component for centralities of 30-80%, 
both of which become large while bubble produc- 
tion becomes smaller as the centrality becomes more 
peripheral. We were able to extend the PBM to 
the PBME and reasonably fit the new quantitative 
RHIC 0-80% y/sWN = 200 Au + Au data[ll for 
the CI and CD correlations. We demonstrated that 
the PBME had a bubble geometry that tracked the 
highest energy density of the different centralities. 
The most central collisions have a symmetry about 
the beam axis. The region of highest energy den- 
sity will also be symmetric about this axis and the 
ring of bubbles centered on this axis is the expected 
geometry nearly identical to that of the PBM. 

As we move to more peripheral collisions the en- 
ergy density decreases. The symmetry of highest 
energy density becomes coupled to elliptic flow and 
is defined by the reaction plane. The region of high- 
est energy density becomes concentrated in the re- 
action plane region. Therefore the ring symmetry 
which applies to central production was broken and 
the bubbles farthest from the reaction plane disap- 
pear. The 12 bubbles for the most central (0-5%) 
reduce to 10 (5-10%), then to 7 (10-20%), and then 
to 5 bubbles (20-30%). Moving to 30-40%, 40-50%, 
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particles from same bubble 




FIG. 8: "(Color online)" The multiplicity(MULT) times the CI correlation for the 0-5% centrality bin as a two 
dimensional vs. A?; perspective plot. This is the part of the CI correlation that resulted from particles emitted 
from the same bubble. 



and 50-60% the bubbles are in the reaction plane 
region with the probability per event of making a 
bubble being 1.5, 0.6, 0.3 (see Table I). 

Both the PBM and the PBME treat the fireball 
in a blast wave model. Bubbles are formed on the 
surface of the blast wave fireball which emit the fi- 
nal state particles at kinetic freezeout. We achieve 
reasonable fits to the quantitative STAR experimen- 
tal analysis of the Charge Independent (CI) and the 
Charge Dependent (CD) correlations within a few 
percent of the total correlations p^. These corre- 
lations considered particles in the pt range 0.8 to 
4.0 GeV/c. The model is also consistent with the 
HBT results in this pt range. The CI correlation 
displays the average structure of the correlated emit- 
ting sources at kinetic freezeout. The CD correla- 
tion has the initial emission correlation of the oppo- 
site sign charge pairs of particles emitted from the 
same space-time region at the time of hadroniza- 
tion. This initial correlation should remain consis- 
tent with the CD correlation of Pythia jets if there 
is no further interactions with other particles af- 
ter hadronization. The analysis above demonstrated 
that at kinetic freezeout the CD correlation in the 
above pt range was consistent with Pythia jets at 
all centralities (see Section VI). Therefore one con- 
cludes that hadronization and kinetic freezeout both 
occur at or very near the surface of the fireball at 
all centralities. Thus the expanding fireball is dense 



and opaque for most centralities at kinetic freeze- 
out. Of course the most peripheral bins due to the 
low material content allow small path length to the 
surface. Both the PBM and its extended present ver- 
sion PBME are surface emission models and thus are 
consistent with this striking experimental feature. 

It is of interest to note that the parton bubble 
model central collision analysis has recently been 
pointed out in [2^ as having features in common 
with glasma flux tubes which evolve from initial 
color glass condensates. This supports the hypothe- 
sis that the bubble substructure can be considered a 
QGP signal and may serve as a key in investigating 
both QGP and glasma effects. 

The persistence of the production of similar sur- 
face bubbles at kinetic freezeout in numbers which 
decrease as the highest energy density decreases as 
the centrality is decreased (going toward the periph- 
eral bins); and their general characteristic of being 
produced in the region of highest energy density im- 
plies the following: 

1) The bubbles represent a significant substruc- 
ture of gluonic hot spots formed on the surface of 
a dense opaque fireball at kinetic freezeout. The 
number of bubbles formed and the energy content 
of each of these substructures is a function of the 
energy density and its extent in space. 

2) Their characteristics, persistence, behavior as 
a function of centrality, and the PBME reasonably 
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FIG. 9: "(Color online)" The multiplicity(MULT) times the CI correlation for the six centralities as six two di- 
mensional A<l) vs. Ar] perspective plots. This is the part of the CI correlation that resulted from particles emitted 
from the same bubble. We see in Fig. 8 and Fig. 9 that the multiplicity scaled correlation grows with centrality. The 
correlation for Ar; remains large, decreases little in the 0-30% centrality bins, and decreases sharply beyond that. 



quantitative fits to the CI and CD data at all cen- 
tralities provide substantial evidence that the bub- 
bles are the final state products of QGP production. 
If sufficiently convincing QGP signatures can be ex- 
tracted from these bubbles then one could eventually 
provide substantial evidence for a QGP. We will be 
investigating this in future work. The bubbles may 



contain relevant information on other topics of in- 
terest (e.g. glasma) which we will also investigate. 
In our future program we plan to utilize the antici- 
pated forthcoming availability of Time of Flight par- 
ticle identification data at STAR to further study 
charged particles correlations with identified parti- 
cles. In the second paragraph from the end of Sec- 
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FIG. 10: The CI correlation bubble signal shape is normalized to one at Ary equals 0.15) vs. Arj for 10° < A^ < 20°. 
Seven centrality bins where there are bubbles present in our model are shown from 50% to 60% increasing to 0% to 
5%. The CI bubble signal Aj; is broad and similar (roughly constant) as a function of Arj for the 0-30% centrality 
range. From Table and text discussion we find that bubbles in this centrality range contains 3-4 partons each and 
that it is the number of partons which determine the Ary width. The 30-60% centrality range show the bubble signal 
A77 width decreases as the number of partons decrease from 3 to 1. 



tion 1.2 of Ref. [Sj we speculated on the possibility 
of applying these model ideas (suitably modified) to 
LHC data when it becomes available. 
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FIG. 11: "(Color online)" The multiplicity(MULT) times the CI correlation for the 0-5% centrality bin as a two 
dimensional Acp vs. Arj perspective plot. This is the part of the CI correlation that resulted from particles emitted 
from all the bubbles. 
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FIG. 12: "(Color online)" The multiplicity(MULT) times the CI correlation for the six centralities plotted as six two 
dimensional A<j!> vs. A77 perspective plots. This is the part of the CI correlation that resulted from particles emitted 
from all the bubbles. We see that the multiplicity scaled CI correlation in Fig. 11 and Fig. 12 from all the bubbles 
behaves similar to that from the CI part which results from a single bubble as shown in Fig. 8 and Fig. 9 with one 
important striking difference. There is a large away sige (> 90°) correlation in the 0-30% centrality range where there 
is large bubble production. This is due to the model conserving momentum between the bubbles and is consistent 
with the experimental observations. 



